Low-temperature heat transport of the layered spin-dimer compound Ba3Mn2 08 
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We report the study on the low-temperature heat transport of Ba3Mn208 single crystal, a layered 
spin-dimer compound exhibiting the magnetic-field-induced magnetic order or the magnon Bose- 
Einstein condensation. The thermal conductivities (k) along both the ab plane and the c axis show 
nearly isotropic dependence on magnetic field, that is, k is strongly suppressed with increasing 
field, particularly at the critical fields of magnetic phase transitions. These results indicate that the 
magnetic excitations play a role of scattering phonons and the scattering effect is enhanced when 
the magnetic field closes the gap in the spin spectrum. In addition, the magnons in the BEC state 
of this materials do not show notable ability of carrying heat. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 
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I. INTRODUCTION 

Low-dimensional or frustrated quantum magnets were 
revealed to exhibit exotic ground states, magnetic exci- 
tations, and quantum phase transitions (QFTs).^' - For a 
particular case of the spin-gapped antiferromagnets, the 
external magnetic field can close the gap in the spectrum, 
which results in a QPT between a low-field disordered 
paramagnetic phase and a high-field long-range ordered 
one. An intriguing finding is that this ordered phase can 
be approximately described as a Bose-Einstein condensa- 
tion (BEC) of magnons.^ Many experimental investiga- 
tions on the critical properties of the BEC-related QPT 
have been carried out, including the characterizations of 
the magnon spectrum, the magnetization, the specific 
heat and the thermal transport, etC/2, In a recent study, 
the heat transport properties of a magnon BEC material, 
NiCl2-4SC(NH2)2 (DTN), were found to display strong 
anomalies at the QPTs and the heat conductivity of the 
BEC state seemed to be much enhanced upon lowering 
temperatures (approaching the absolute zero).** This re- 
sult shows an analogy between the magnon BEC and the 
superfiuid of ''He in the aspect of the ability of transport- 
ing heat. However, one notable facet is that magnons act 
as heat carriers only in the direction of the spin chains of 
this compound, whereas they only scatter phonons in the 
transverse directioni^ So the exchange anisotropy may 
play the key role in the heat transport of magnetic exci- 
tations. A later experimental work confirmed the main 
features of transport properties of DTN, but an alter- 
native picture based on the mass renormalization and 
impurity scattering effects was proposed to explain the 
thermal transport data.^ In any case, the QPTs associ- 
ated with the magnon BEC are believed to significantly 
affect the heat transport properties and the low-energy 
magnetic excitations provide a substantial contribution 
to the heat transport. To get the general principals of 
the heat transport in the magnon BEC state, we need 



to carry out systematic studies on more members of the 
magnon BEC materials. 

Ba3Mn208 (BMO) is an 5 = 1 quantum spin-dimer 
system exhibiting the BEC of spin degrees of freedom. It 
crystallizes in the rhombohedral i?3m structure with the 
pairs of Mn^+ 3cf ions (5 = 1) arranged vertically on the 
hexagonal layers.^ The spins of Mn^+ ions of each pair 
are coupled antiferromagnetically to form spin dimers.— ^— 
Neutron scattering results indicated an intradimer ex- 
change energy Jq = 1.642 meV, an interdimer coupling 
between Mn ions in the same plane J2 — J3 = 0.1136 meV, 
an interdimer coupling between Mn ions in the adjacent 
planes Ji = —0.118 meV, and the next nearest neigh- 
bor interdimer interactions between bilayers (along the c 
axis) J4 = —0.037 meVM-^H The strong in-plane spin- 
dimer interaction results in a spin-singlet ground state, 
with a spin gap of 1.05 meV to the lowest triplet state 
and a second larger gap to the quintuplet state. BMO 
displays a peculiar phase diagram in magnetic field when 
the Zeeman effect lowers the Sz — I triplet states and 
the Sz = "2 quintuplet statesJ^ Corresponding to the 
closures of two spin gaps with increasing magnetic field, 
there are two sequential magnetically ordered states and 
four quantum critical fields: Hd = 8.7 T, Hc2 = 26.5 T, 
Hc3 = 32.5 T, Hci = 47.9 T for H \\ cJ^^^^ The pecu- 
liarity is that the first ordered phase for Hd < H < Hc2 
and the second one for Hcs < H < Hd are described 
as the triplet condensation and the quintuplet conden- 
sation, respectively." Furthermore, BMO was found to 
have a weak single-ion uniaxial anisotropy, characterized 
by \D\ = 0.032 meV.— The phase diagram for _ff _L c is 
a bit more complicated, with two phases of the triplet 
condensationFi^"— It was revealed that both the single- 
ion anisotropy and the geometry frustration play crucial 
roles in determining the phase diagram.— 

In this work, we study the thermal conductivity (k) of 
BMO single crystal to probe the field-induced QPTs and 
the role of magnetic excitations in the heat transport. It 
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FIG. 1: (Color online) Low-temperature specific heat of 
Ba3Mn208 single crystal in zero field. The solid line indi- 
cates the fit of data using formula ([T}. Inset: the data in a 
broader temperature range from 0.4 to 20 K. 



is found that both Kab and strongly suppressed in 

magnetic fields, particularly at the field-induced QPTs, 
which demonstrates that at very low temperatures the 
magnons mainly act as phonon scatterers even in the 
long-range ordered state. The present results indicate 
that magnon-BEC state may not necessarily exhibit large 
thermal conductivity and the ability of magnetic excita- 
tions transporting heat is determined by the magnetic 
structure and the anisotropic spin exchange. 



II. EXPERIMENTS 



High-quality single crystals of Ba3Mn2 08 are grown by 
a slowly cooling method using NaOH as a solvent The 
crystals have a hexagonal-platelet-like shape with large 
sizes up to 5 x 5 x 2.5 mm"^. The largest surface is the ab 
plane and the thickness is along the c axis, determined by 
using the x-ray diffraction and the Laue photograph. The 
specific heat is measured by the relaxation method in the 
temperature range from 0.4 to 20 K using a commercial 
physical property measurement system (PPMS, Quan- 
tum Design). Both the a6-plane and the c-axis thermal 
conductivities (nab and Kc) are measured using the con- 
ventional steady-state technique. Particularly, the data 
at low temperature and in high magnetic field are taken 
by using a "one heater, two thermometer" method in 
a ■^He refrigerator and a 14 T magnet Note that 

we can only probe the QPTs at the first critical field of 
BMO, due to the limitation of the available field. 



III. RESULTS AND DISCUSSION 

It is known that, in a magnetic insulator, heat can be 
carried by phonons and magnetic excitations, and the 
interaction between them usually induces scattering on 
heat carriers, which is a negative effect on the heat trans- 
port. To analysis the heat transport properties quantita- 
tively, it is necessary to know the phononic specific heat 
from the experiments. Figure 1 shows the low-T specific 
heat of the Ba3Mn208 single crystal. At ~ 6 K, there ap- 
pears a broad shoulder-like feature, which is apparently 
of magnetic origin. Below 5 K, the specific heat decreases 
quickly with temperature and it does not show a de- 
pendence of the phonon specific heat. The same result 
has been reported in an earlier worki^ The low-T specific 
heat was found to be able to described by a formula for 
the gapped magnetic excitations?^ 



C = nR{A/Tfe- 



-A/T 



(1) 



where h is the number of excited states per spin dimer, 
R is the gas constant, and A is the energy gap of mag- 
netic excitations. The low-T data can be fitted rather 
well to Eq. ([T]) with parameters n — 1.49 and A = 14.2 
K, except for the very low-T data (< 1 K)4 Note that 
the size of gap is quite consistent with those from other 
measurements jiiii However, it is difficult to separate the 
phononic specific heat because of the significant contri- 
bution from the magnetic excitations. Actually, adding a 
T'^-term to Eq. (H]) cannot achieve perfect fitting to the 
data, particulary for the subKelvin data. 

Figure 2 shows the temperature dependencies of Kab 
and Kc of BMO single crystals in zero and 14 T fields, 
which are applied along either the ab plane or the c 
axis. Apparently, the zero-field thermal conductivities 
show rather weak anisotropy at high temperatures. Each 
k{T) curve exhibits a phonon peak at 18 K with the peak 
values nearly the same for Kab and Kc- The T^^ depen- 
dence of k{T) at high temperatures is the characteristic of 
phonon heat transport dominated by the phonon-phonon 
Umklapp scatteringjii At sub-Kelvin temperatures, both 
Kab{T) and Kc{T) show an approximate T^'^ dependence, 
which indicates that the boundary scattering limit of 
phonon thermal conductivity is approached or that the 
microscopic scattering on phonons is negligibleJ^ Be- 
sides these usual behaviors belonging to the phonon heat 
transport, there is a "shoulderlike" feature in the zero- 
field k{T) at 6.5 K, which is likely due to some kind of 
resonant scattering on phonons>ii Applying 14 T mag- 
netic field leads to a strong suppression of low-T ther- 
mal conductivity and the disappearance of the "shoul- 
derlike" feature. Therefore, the phonon resonant scat- 
tering in zero field must be of magnetic origin. Similar 
result has been found in another spin-gapped material, 
DTNji On the other hand, there is no evidence for the 
magnon heat transport in BMO, which is reasonable at 
very low temperatures considering the negligible magnon 
excitations due to the finite spin gap. In high magnetic 
field, the strong suppression of phonon heat transport 
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FIG. 2: (Color online) Temperature dependencies of the ther- 
mal conductivities Kab and Kc of Ba3Mn208 single crystals in 
zero and 14 T fields. The dashed lines indicate the approxi- 
mate T^'^ and dependencies of the thermal conductivity 
at very low temperatures and at high temperatures, respec- 
tively. 



in a broad temperature range is apparently due to the 
enhancement of magnetic scattering, which can also be 
seen from the following magnetic-field dependencies of k. 
Another peculiar feature of 14 T k(T) data is the small 
jumps of conductivity at low temperatures, which are 
related to the field-induced antiferromagnetic (AF) tran- 
sitions. Note that even with these increases of k, at the 
phase transitions, the lower-T thermal conductivities are 
still much smaller than the zero-field values, indicating 
that the phonons are still strongly scattered by magnetic 
excitations in the field-induced AF state. 

Figure 3 shows the detailed magnetic field dependen- 
cies of Kc and Kah for both _ff || a6 and 77 || c. It is 
clear that they are nearly isotropic on the field direction. 
The overall behavior of K(-ff) is that the magnetic field 
strongly suppresses thermal conductivity, but at very 
low temperatures, i.e., at 0.36 and 0.5 K, k shows up- 
turn above some critical fields. Since the zero-field ther- 
mal conductivity of BMO is purely phononic, the field- 
induced suppression of k is clearly a result of scattering 
effect on phonons by some magnetic excitations. In this 
regard, the phonon scattering by paramagnetic moments, 
which are commonly existing as impurities or spin vacan- 
cies, is one of the unavoidable effectSfi2r2S However, the 
field dependencies of k shown in Fig. 3 are rather dif- 
ferent from the standard one of paramagnetic moments 
scattering.—"^ The main origin for the behaviors is 
therefore the phonon scattering by magnon excitations. 
The low-field quantum disordered phase has a finite spin 
gap, which can be weakened by the Zeeman effect. At 
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FIG. 3: (Color online) Magnetic-field dependencies of ther- 
mal conductivity of Ba3Mn208 single crystals at low temper- 
atures. 



a fixed temperature, with increasing the magnetic field, 
the spin gap decreases and the number of low-energy 
magnons increases quickly, which can strongly scatter 
phonons and leads to a significant suppression of k. 

A notable feature of isotherms is the upturn of 

K in high magnetic field and at very low temperatures. 
This increase of k is obviously related to the field-induced 
AF order since the transition fields are found to be co- 
incided with the phase boundary of high-field ordered 
phase, i.e., the magnon BEC state. There are two possi- 
ble reasons for the increase of k in the BEC state. First, 
the phonon scattering by magnons is weakened when the 
magnon BEC occurs. Second, the magnons in the BEC 
state have a positive effect on n, by transporting the heat 
directly. If the magnons transport the energy, due to the 
anisotropy of the exchange and the low-dimensionality of 
the spin structure, the field dependence of k should ex- 
hibit obvious anisotropy on the direction of heat current 
with the same fieldi^"— However, it can be seen from 
Fig. 3 that the upturn behaviors actually show a weak 
anisotropy on the direction of magnetic field rather than 
on that of the heat current. This means that the upturn 
of thermal conductivity in the BEC state is a result of 
the weakening of magnon-phonon scattering. In other 
words, although the magnon scattering on phonons are 
generally strong at high fields, it is more significant at 
the critical fields because of the strong spin fiuctuations. 

Although the above discussions are quite clear in catch- 
ing the main physics of the observed heat transport prop- 
erties, it is useful to try a more quantitative analysis on 
the experimental data. First, since the phonons are the 
only one type of heat carriers, the zero-field K(r) data are 
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likely to be fitted by a classical Debye model of phonon 
thermal conductivity^^'^"^ 



Qu/T 



1 



,T{u,T)dx, (2) 



in which uj is the phonon frequency, x — huj/ksT is di- 
mensionless, and 1/t(cj, T) is the phonon relaxation rate. 
The relaxation time is determined by 

= v/L + Alj^ + BTuj^ expi-Oo /bT) + r"^, (3) 

which represent the phonon scattering by the grain 
boundary, the point defects, the phonon-phonon 
Umklapp scattering, and the resonant scattering, 
respectivelyfii The average sound velocity v can be cal- 
culated from the Debye temperature using the formula 
©D — w(?i/fcs)(67r^n)^/'^, where n is the number density 
of atoms. For BMO, however, the Debye temperature 
cannot be simply obtained from the phonon specific-heat 
data, as the data in Fig. 1 indicate. For an estimation, we 
refer the specific- heat data of Ba3V208r^ a nonmagnetic 
material having the same crystal structure as BMO. The 
phonon specific-heat coefficient /3 — 6.06 x 10^* J/K*mol 
can be easily obtained from a T'^-fitting to the low-T 
datai^ With this /3 value, the Debye temperature and 
the sound velocity of BMO are calculated to be 346 K 
and 2900 m/s, respectively. The parameter L describing 
the boundary scattering is the the averaged sample width 
of the samples. The other parameters A, B and b are free 
ones. 

One difficulty of getting a precise calculation using Eq. 
^ is the description of t^^. In BMO, the magnetic 
excitations can strongly scatter phonons. Considering 
the energy dispersion of magnons, it is very difficult to 
get a transparent formula for r^^. For a simple analysis, 
one can neglect the magnon dispersion and assume that 
the resonant scattering of phonons occurs via a singlet- 
triplet excitations of dimerized states. Thus, the resonant 
scattering rate can be expressed as 



C 



FiT), 



(4) 



where C is a free parameter while F{T) describes the dif- 
ference of thermal populations of the excited triplet and 
the ground singlet states (the contribution of a higher- 
level excitations of quintuplet states is negligibly small). 
It is known thai^^ 



F{T) = 1 - 



1 - exp{-A/T) 
l + 3exp{-A/Ty 



(5) 



with A hujo the energy gap of magnetic spectrum, 
which is about 14.2 K from the specific-heat data. 

The zero-field thermal conductivity data are able to be 
fitted by using the above formula very well. As an ex- 
ample, the zero-field data from Fig. 2(a) and the fitting 
curve are displayed in Fig. 4(a). The best fitting param- 
eters are A ^ 1.3 x lO"^"^ s^, B = 3.0 x IQ-^o K^^s^, b 
= 6.0, and C = 1.6 X 10^ s^^ 
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FIG. 4: (Color online) Comparison of the Kab{T) data (a) 
and the Kab{H) data (b) with the calculation results using 
the Debye model. 



Second, the magnetic-field dependence of k can also be 
calculated using the Debye model. The effect of magnetic 
field is to induce the Zeeman splitting of the degenerate 
triplet states and a linear decrease of the energy gap with 
field, that is. 



AiH)^A-gfiBH/kB. 



(6) 



The Lande factor was known to be ^ 2.00.^^'^^ A typical 
comparison between the k{H) data and the calculations 
is shown in Fig. 4(b). Although there is a qualitative 
agreement between the experimental data and the calcu- 
lation, the quantitative difference is quite large. In par- 
ticular, at low fields the experimental data show much 
stronger suppression than the calculated result. It is ex- 
pectable from the classical Boltzmann distribution that 
the magnon excitations are negligibly weak at subKelvin 
temperatures when the spin gap is larger than 10 K. The 
experimental data suggest that the magnon excitations, 
which may be because of the quantum fluctuations, are 
much stronger than the thermal excitation produces. An- 
other reason for the significant discrepancy is likely re- 
lated to the simple assumption of the resonant phonon 
scattering by the two-level magnetic excitations. Actu- 
ally, the magnon dispersions of BMO are not very weak, 
which makes the scattering between phonons and mag- 
netic excitations inelastic and much more complicated. 
A precise calculation based on the magnon-phonon scat- 
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FIG. 5: (Color online) Temperature dependencies of the ther- 
mal conductivity measured in magnetic field above the critical 
field Hci ■ The arrows indicate the critical temperatures where 
ft starts to increase. 



tering calls for the details of both the phonon spectra and 
the magnon dispersions .^^^^ 

It is useful to compare these results with those of DTN, 
in which the magnetic field along the c axis (the spin- 
chain direction) can close the spin gap and drive the AF 
phase transition. In DTN, the heat transport along the c 
axis or perpendicular to it present remarkably anisotropic 
behaviors with increasing field, that is, Kc show sharp 
peaks at Hd while Kab show dips at the same position.— 
Those demonstrated that magnons act as the heat car- 
riers along the c axis and the phonon scatterers along 
the ab plane, respectively."* Furthermore, the Kc of DTN 
in the AF state (for Hd < H < Hc2) keep increasing 
with lowering the temperature, also pointing to a poten- 
tially large heat transport in the magnon BEC phase. ''■^ 
The present k{H) behaviors of BMO for both Kab and 
Kc are very similar to those of tiab{H) in DTN, which 
suggests that the magnetic scattering on phonons are en- 
hanced upon approaching the critical fields. In contrast 
to the case of DTN, the magnetic excitations of BMO do 
not contribute to transporting heat substantially along 
either the ab plane or the c axis. Therefore, whether the 
magnon BEC state of the spin-gapped compounds can 
have strong ability of transporting heat is mainly depen- 
dent on such factors as the nature of the magnetic struc- 
ture and the exchange anisotropy etc. The field-induced 
AF state or magnon BEC is not the sufficient condition 
for observable magnon transport. 

The above are, however, not the whole story because 
there is only one transition in each k{H) isotherm for 
both II c and H \\ ab, whereas the former experiments 
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FIG. 6: (Color online) The transition points of k{T) curves 
from Fig. 3. Temperature-field phase diagram, obtained by 
heat capacity and magnetocaloric effect measurements from 
Ref. IT2 . is shown for comparison. I and II indicate two dif- 
ferent phases of the triplet condensate.— 



revealed a single magnetic transition and two distinct 
transitions, respectively. For more precise description, we 
show in Fig. 5 a series of k{T) curves for magnetic fields 
varying between 10 and 14 T. It is found that below 10 T, 
there is no any anomaly in the k{T) curves down to 0.3 K. 
Above 10 T, a jump in k{T) shows up and becomes bigger 
with increasing field up to 14 T. From these data, we can 
obtain the transition points on the H — T phase diagram 
and compare them with the boundaries of the magneti- 
cally ordered phases determined by heat capacity, mag- 
netocaloric effect, and cantilever torque measurements 
as shown in Fig. 6. It is interesting that the jump-like 
anomalies of k occur at the boundary of phase I, that is, 
the thermal conductivity increases at the phase transi- 
tions to phase I from either the low-field disordered state 
or the antiferromagnetically ordered phase II. This phe- 
nomenon has a good correspondence to the specific-heat 
data,*^ which shows a sharp lambdalike peak and a much 
weaker and less divergent peak at the phase transitions 
from the disordered state or phase II to phase I and the 
disordered state to phase II, respectively. It is likely that 
the phase transition associated with lambdalike specific- 
heat anomaly has more significant critical fluctuations 
that strongly scatter phonons. 
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IV. SUMMARY 

We study the heat transport of Ba3Mn208 single crys- 
tals at very low temperatures and in high magnetic fields 
to probe the field-induced magnetic phase transitions 
and the role of magnons in the transport properties. In 
zero field, the low-T thermal conductivity shows a purely 
phononic transport behavior. In applied strong field, the 
spin gap is gradually diminished and the magnetic excita- 
tions are populated, which induces an enhanced magnetic 
scattering on phonons and an overall strong suppression 
of thermal conductivity with increasing field. Moreover, 
at very low temperatures where the closure of spin gap 



results in a long-range ordered AF state, the phonons 
are more significantly scattered at the phase transitions. 
There is no evidence showing the sizeable magnetic heat 
transport in this magnon BEC compound. 
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